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Abstract 

We study the influence of nonlocality in the interaction on two spin one pairing 
patterns of two-flavor quark matter: the anisotropic blue color paring besides the 
usual two color superconducting matter (2SCb), in which red and green colors are 
paired, and the color spin locking phase (CSL). The effect of nonlocality on the gaps 
is rather large and the pairings exhibit a strong dependence on the form factor of 
the interaction, especially in the low density region. The application of these small 
spin-one condensates for compact stars is analyzed: the early onset of quark matter 
in the nonlocal models may help to stabilize hybrid star configurations. While the 
anisotropic blue quark pairing does not survive a big asymmetry in flavor space as 
imposed by the charge neutrality condition, the CSL phase as a flavor independent 
pairing can be realized as neutral matter in compact star cores. However, smooth 
form factors and the missmatch between the flavor chemical potential in neutral 
matter make the effective gaps of the order of magnitude ~ 10 keV, and a more 
systematic analysis is needed to decide whether such small gaps could be consistent 
with the cooling phenomenology. 
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1 Introduction 


The investigation of color superconducting phases in cold dense quark matter 
has attracted much interest in the last years since those phases could be rel¬ 
evant for the physics of compact star cores [1]. Nevertheless, it became clear 
that large spin-0 condensates, like the usual two-flavor superconductor (2SC), 
although having large pairing gaps (~ 100 MeV) and therefore a direct influ¬ 
ence on the equation of state (EoS), may be disfavored by the charge neutrality 
condition if not unusually strong diquark coupling constants are considered 
[2,3,4], 

Alternatively, spin-1 condensates are being investigated [5,6,7], Due to their 
smallness, their influence on the EoS is negligible but they could strongly af¬ 
fect the transport properties in quark matter and therefore have important 
consequences on observable phenomena like compact star cooling, see [8], A 
recent investigation of neutrino emission and cooling rates of spin-1 color su¬ 
perconductors constructed for conserved total angular momentum allows for 
color-spin locking, planar, polar, and A phases [9], However, none of these 
phases fulhlls the requirements of cooling phenomenology that no ungapped 
quark mode should occur on which the direct Urea process could operate and 
lead to very fast cooling in disagreement with modern observational data [12], 

In the present work, we consider spin-1 pairing patterns different from the 
above mentioned, like the anisotropic third color (say blue) quark pairing 
besides the usual 2SC phase (2SCb) [10] and the s-wave color-spin locking 
phase (CSL) [11], which have been introduced within the NJL model whereby 
small gaps in the region of some fractions of MeV have been obtained. Such 
small gaps could help to suppress efficiently the direct Urea process in quark 
matter and thus to control the otherwise too rapid cooling of hybrid stars 
[12,13], 

One important feature is that the form of the regularization, i,e, via a sharp 
cut-off or a form factor function, is expected to have a strong impact on the 
resulting gaps due to the sensitive momentum dependence of the integrand in 
the gap equation [10], Especially the behavior of quark matter in the density 
region of the suspected deconhnement transition plays a crucial role for deter¬ 
mining the stability of compact star conhgurations. Models with a late onset 
of quark matter could eventually lead to unstable hybrid star conhgurations. 

For example, in [11] it has been shown within the local NJL model how a pos¬ 
sible pairing pattern for compact star matter could be described which fulhlls 
the constraints from compact star cooling phenomenology. These require that 
all quark species be paired and the smallest pairing gap be of the order of 
10 keV to 1 MeV with a decreasing density dependence, A caveat of the NJL 
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model quark matter is, however, that a stable hybrid star can be realized only 
marginally, see [14], 

The advantage of nonlocal models is that they can describe the regularization 
of the quark interaction via form factor functions and therefore represent it 
in a smoother way, especially for low densities [15]. For the case of the 2SC 
phase it has already been shown that the effect of nonlocality in the low density 
region is rather large and the pairing exhibits a strong dependence with the 
form factor of the interaction [2,16]. Moreover, the early onset of quark matter 
for the dynamical chiral quark model, in contrast to the NJL model, might 
help to stabilize hybrid star configurations. 

As it has been shown in [17] within a nonlocal generalization of the NJL model 
[15], stable hybrid stars with large qnark matter cores can be obtained. In order 
to describe the properties of these stars consistently, inclnding their cooling 
phenomenology, the description of diqnark pairing gaps as well as emissivities 
and transport properties for the above motivated CSL phase should be given 
also within a nonlocal quark model. 

As a hrst step in this direction we will provide in the present paper the spin-1 
pairing gaps for a nonlocal, instantaneous chiral quark model under nentron 
star constraints for later application in the cooling phenomenology. We give 
here a hrst discussion of the inhuence of nonlocality of the interaction in 
momentnm space when compared to the NJL model case and discuss the 
possible role of spin-one condensates for compact star applications. 

The paper is organized as follows: in Section 2 we briehy review how the NJL 
model is modihed when nonlocality is introdnced using a three-dimensional 
nonlocal chiral qnark model; in Sections 3 and 4 we present the nonlocal 
version of the anisotropic blue color paring (2SCb) besides the usnal two- 
havor color snpercondncting (2SC) phase and of the color-spin locking phase 
(CSL), respectively. In Section 5 we present preliminary results for neutral 
matter in compact stars for the Gaussian form factor and discuss whether 
requirements for hybrid star cooling phenomenology conld be met. Finally, in 
Section 6, we draw the Conclusions. 


2 Nonlocal chiral quark model 


We investigate a nonlocal chiral model for two-havor quark matter in which 
the qnark interaction is represented in a separable way by introducing form 
factor functions g{p) in the bilinears of the cnrrent-current interaction terms 
in the Lagrangian [2,15,16,17]. It is assnmed that this four-fermion interaction 
is instantaneous and therefore the form factors depend only on the modulus 
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of the three momentum p = |p|. 

In the mean held approximation the thermodynamical potential can be eval¬ 
uated and is given by 


= J-^;^^TT\n{^S \iuJn,p)) + V, ( 1 ) 

where the sum is over fermionic Matsubara frequencies Un = {2n+l)7rT and V 
is the quadratic contribution of the condensates considered. The specihc form 
of V in dependence on the order parameters 0 for chiral symmetry break¬ 
ing and A for color superconductivity in the corresponding diquark pairing 
channels will be given below in Section 3. 


In our nonlocal extension the inverse of the fermion propagator in Nambu- 
Gorkov space is modihed in comparison to the NJL model case by momentum 
dependent form factor functions g{p) as follows, 


S-\p) 


^ ^ + /i7° — M{p) g{p)^ 

^ -^(p)A^ pi- - M{p) 


( 2 ) 


where fi is the chemical potential matrix and the elements of M{p) = diag{Mj(p)} 
are the dynamical masses of the quarks given by 


Mf{p) =mf + g{p)(l)f . (3) 

The matrix A represents the order parameters for diquark pairing which will 
be made explicit in Sects. 3 and 4, respectively. 

In (2) and (3) we have introduced the same form factors g{p) to represent the 
nonlocality of the interaction in the meson [qq) and diquark [qq) channels. In 
our calculations we use the Gaussian (G), Lorentzian (L) and cutoff (NJL) 
form factors dehned as 


gG{p)=exp{-p^/Al) , (4) 

gM = [1 + ( 5 ) 

5'njl(p) = ^(1 — p/Anjl) • (6) 

The parameters for the above form factor models used in this work are pre¬ 
sented in Tab. 1. They have been hxed by the pion mass, pion decay constant 
and the constituent quark mass M = M(0) at T = /x = 0. In order to estimate 
the effect of the nonlocality on the results relative to the NJL model, we used 
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go, Ql and Qnjl with parameters fixed such that M = 380 MeV, see Tab. 1. 
For details of the parameterization, see [22], 


Form Factor 

Notation 

A [MeV] 

G A2 

m[MeV] 

Gaussian 

9G 

786.7 

4.12 

2.50 

Lorentzian 

9L 

637.2 

2.76 

2.59 

NJL 

9njl 

596.1 

2.36 

5.54 


Table 1 

Parameter sets for the nonlocal chiral quark model {go, Ql ) and for the NJL model 
(gNJi)', for all M(0) = 380 MeV is fixed. 


The stationary points of the thermodynamical potential (1) are found from 
the condition of a vanishing variation 

= 0 (7) 

with respect to variations of the order parameters. Eq. (7) defines a set of gap 
equations. Among the solutions of these equations the thermodynamically 
stable state corresponds to the set of order parameter values for which hi has 
an absolute minimum. 


3 Anisotropic blue quark pairing for nonlocal model 


First we consider the 2SCb phase in which two of the three colors (e.g. red 
r and green g) pair in the standard spin-0 isospin singlet condensate (2SC 
phase) and the residual third color (consequently, it is blue b) pairs in a spin -1 
condensate (symmetric in Dirac space, symmetric in color, antisymmetric in 
flavor) [10]. The matrix A in the inverse quark propagator ( 2 ) for the 2SCb 
phase is then given by 

= A{'y,T2X2){Sc,r + 6,,g) + A'(a°3 r2 Pt%,b , (8) 

where T 2 is an antisymmetric Pauli matrix in the flavor space and A 2 an an¬ 
tisymmetric Gell-Mann matrix in the color space; = |[ 7 °, 7 ^] and = 
A — ^ is the projector on the third color. If A' 7 ^ 0, this can be understood 
as a nonzero third component of a vector in color space which breaks the 
0{3) rotational symmetry spontaneously. The blue quark pairing is therefore 
anisotropic. We consider first symmetric matter and thus we take the quark 
chemical potentials as Pu = Pd = p. 
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The thermodynamical potential is given by 




3 

-4E 

i=l ' 


4Gi 

d?p 

(27r)3 


4G2 I6G3 
e- + e: 


^ + T In (1 + /^) + T In (1 + 6“^^ /^) 


( 9 ) 


where the coupling constants Gi,G 2 ,G 3 follow the relation given by the in- 
stanton induced interaction [18] 


Gi : ^2 : G 3 = 1 : 3/4 : 3/16 . (10) 

The dispersion law for the paired quarks (r, g) is given by 

EUp) = E^(p) = ^(€ T + 9"(p)l AP , (11) 

where e = is the free particle dispersion relation and M = = 

Md. 

For the anisotropic pairing of the blue quarks the dispersion relation can be 
written as 

= V(^efrThefr)^ + c/^(p)|A(,g|2 , (12) 

where the effective variables depend on the angle 9, with cosO = ps/lp], and 
are dehned as 


el^ = p^ + Ml^, 

(13) 


(14) 

/^eff 


Fefr = F^ + /(p)|Af sin^^. 

(15) 

A;g|2 = |Af(cos20 + ^sin2 0). 

(16) 


Feff 


The dispersion relation Ef is, therefore, an anisotropic function of p and 
therefore the calculation of (9) should be performed as an integral over the 
modulus of IpI and over the angle 9. Ef has a minimum if 9 = 7r/2 and 
vanishes if M = 0 or A' = 0. 

As it has been pointed out in [10], the gap equations for A and A' are only 
indirectly coupled by their dependence on M. Therefore, since the equation 
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for A' nearly decouples if M is small, we can illustrate the anisotropic con¬ 
tributions to the thermodynamical potential hi fixing the variables (/i, T) and 
the order parameters (0, A) and varying the angle 6. For this purpose, we 
consider 

^ d{cos |cos6» , (17) 

and in Fig. 1 we plot I^ose as a function of the gap A'. As 9 increases 

from 0 to 7r/2 the position of the minimum of moves to lower values 

of the gaps A'. The value of A' that minimizes the thermodynamical potential 
is found once the integration over the angle 6 is performed. 



Fig. 1. Anisotropic contributions |cose to the thermodynamical potential 

^ 2 SCb g fuHction of tFc spin-1 gap A' for fixed (p = 500 MeV, T=0) and (0 = 20 
MeV, A = 106 MeV). As 9 increases the position of the minimum of |cos6» 

moves to lower values of the gaps Ah The low gap region is zoomed in the inset 
figure on the bottom left. For (solid line) we obtain that the minimum is 

placed at A' ~ 2 MeV/fm^. The NJL parameterization is considered. 


3.1 Gap equation solutions 


We search for the stationary points of (7) respect to the order parameters 
solving the gap equations 
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5f]2SCb 5f]2SCb 


(50 


(5A 


6A’ 


= 0 


(18) 


using the dispersion relations (11) and (12). The results that are shown in this 
work are obtained for T = 0. 


In Fig. 2 we show the chiral gap 0, the 2SC diquark gap A and the spin-1 
pairing gap of the blue quarks A' as functions of the quark chemical potential 
fi. The gaps A' are strongly density-dependent rising functions and typically 
of the order of magnitude of keV, e.g., for a hxed /i at least two orders of 
magnitude smaller than the corresponding 2SC gaps. These small gaps are 
very sensitive to the form of the regularization and to the parameterization 
used. Obviously, also the onset and the slope of the superconducting phases 
depend strongly on the parameters used. 



Fig. 2. Chiral gap 0, 2SC diquark gap A (upper panel) and spin-1 pairing gap of 
the blue quarks A' (lower panel) as a function of the quark chemical potential p for 
the NJL model and for the Gaussian form factor. 

The effect of the nonlocality on the results is shown in the lower panel of Fig. 2: 
the smoothness of the Gaussian form factor reduces the A' gaps dramatically. 












For this case, the blue quark pairing gaps are about two orders of magnitude 
lower than the corresponding NJL model (both with hxed M = 380 MeV). 

The dependence of the results on the parameters used is also rather strong and 
nonlinear as it is shown in Fig. 3 using, as example, the Gaussian form factor. 
When the coupling constant G 3 is doubled (dash-dotted line) the resulting 
gaps increase between two and three orders of magnitude, depending on the 
chemical potential. 



Fig. 3. Effect of the coupling constant G 3 on the spin-1 blue quarks energy gaps A' 
for the Gaussian form factor. The dash-dotted line corresponds to the case when 
the coupling constant G 3 is doubled. 


Since these small gaps are practically negligible in comparison to usual 2SC 
gaps, they would have no influence on the equation of state. On the other hand, 
it is well known that even small pairing energies could play an important role in 
the calculation of transport properties of quark matter for temperatures below 
the order of magnitude of the gap parameter. Nevertheless, it is unlikely that 
the blue quark pairing could survive the compact star constraints: the charge 
neutrality requires that the Fermi seas of the up and down quarks should differ 
by about 50-100 MeV and this is much larger than the gaps that we obtain 
for these condensates (~ keV) in the symmetric case. 

In the following section we present the nonlocality effects on a flavor symmetric 
spin -1 pairing channel being a good candidate to survive the large mismatch 
between up and down quark Fermi seas in charge neutral quark matter. 
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4 Color-spin locking (CSL) phase for a nonlocal chiral model 


In the s-wave CSL phase introduced in Ref. [11], which differs from the CSL 
phase in [7], each condensate is a component of the antisymmetric anti-triplet 
in the color space and is locked with a vector component in the spin space. 
In the present paper, we study a nonlocal generalization of the CSL pairing 
pattern of Ref. [11] and consider the matrix A in (2) for the CSL channel as 

= A/(73A2 -f 71-^7 + 72 A 5 ) • (19) 

The thermodynamical potential can be decomposed into single-flavor compo¬ 
nents 

{f,.,}) = ^ (20) 


where the contribution of each flavor is 


nfyr.. I- ^ 11 ^ 2T1„ (1 . e- V./ 7 , 


( 21 ) 


The ratio of the two coupling constants 


Gi : i/. = 1 4 . ( 22 ) 

O 

is obtained via Fierz transformations of the color-currents for a one-gluon 
exchange interaction. To derive the dispersion relations Ef^k we follow [11] 
and we extend the expressions for the nonlocal model introducing the form 
factors to modify the quark interaction. We obtain that Ef.^ could be brought 
in the standard form 


- (^/,eff T + |^5'^(p) (23) 

if the effective variables are now dehned as: 

+ ^/,efr 1 
Afyeff = Mf{p) ) 

h/,eff 

/^7eff = /^? + |A/lV(p), 
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(24) 

(25) 

(26) 



(27) 


A/,efr - 


Mfjp) 


A, 


For Ef-k, k = 3...6, the particle and the antiparticle branches split 


-^/;3,5 - (^/ “ h/) +«/;3,5|A/| Q (p) , 

EfA.e = (<?/ + + 0/i4,6lA/p9'^(p) , 


(28) 

(29) 


where the momentum-dependent coefficients af-k, k = 3...6 are given by 


«/;3,5 - 


*^/;4,6 


5 P" 

± A 

(l ' 

]\s^Ep) 


\ 



5+ ^ 

± A 

(l+ ^ 

f + sWrt 



V ^/h/y 



(30) 

(31) 


and 

e)=p^ + M]{p) . 

We solve the gap equations 


(32) 




f 


= 0 


(33) 


( 50 / 6Af 

and present the results of the global minimum of in the next subsection. 


4-1 Gap equation solutions for each flavor 


The mass gaps and the CSL gaps are shown in Fig. 4 for different form factors 
of the quark interaction. 

The CSL gaps are strongly /i/-dependent rising functions in the domain that is 
relevant for compact star applications. There is a systematic reduction of the 
CSL gaps as the form factors become smoother (from NJL to Gaussian) and 
the condensates in the nonlocal extension are at least one order of magnitude 
smaller than in the NJL case. 
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325 350 375 400 425 450 475 500 

[MeV] 


Fig. 4. The dependence of the chiral gap (j) and the CSL pairing gap Aj on the 
chemical potential /i/ for different form factors of the nonlocal interaction and for 
the NJL model. 

Note that especially the low-density region is qnalitatively determined by the 
form of the interaction. Since the nonlocality also affects the chiral gap, the 
breakdown of which is a prereqnisite for the occnrrence of color snpercon- 
dncting phases, we observe for the nonlocal models an earlier onset of the 
snpercondncting qnark matter (/^/,crit < 350 MeV) than in the NJL model 
cases. As it has been pointed ont in [17], the position of the onset is crncial to 
stabilize hybrid star conhgnrations. In general, NJL models present a later on¬ 
set than nonlocal ones [16] and might disfavor the occurrence of stable hybrid 
star configurations with a quark matter core [19]. 

In Fig. 5 we plot the effective CSL gaps setting the explicit dependence of the 
form factor g{p) = 1 in (27) in order to compare the order of magnitude of 
them with sharp cut-off models. We obtain that the Gaussian is an increasing 
function of the chemical potential from approximately 15 keV to 35 keV and 
the Lorentzian is nearly constant of the order of 80-90 keV. Both exhibit gaps 
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that are much smaller than the corresponding NJL ones which are in the range 
Aj,efr ^ 300-200 keV. 



Fig. 5. The effective CSL pairing gap ^f,eS for different form factors as a function 
of the chemical potential pj. 


5 Results for matter in compact stars 


Since the CSL pairing is symmetric in flavor, we can easily construct electri¬ 
cally neutral quark matter in /3-equilibrium for compact star applications. We 
consider stellar matter in the quark core of compact stars consisting of {u, d} 
quarks and leptons {e, Ue, he, /i, The particle densities rij are conjugate 

to the corresponding chemical potentials fij according to 


on 


Hj = 


d^j 


00 AoT 


(34) 


where the index j denotes the particle species. We consider matter in (3- 
equilibrium with only electrons and since we assume that neutrinos leave the 
star without being trapped (/i^^ = —/ii/^ = 0, /ip^ = —/Xi/^ = 0) 
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/^e — • 


( 35 ) 


We are interested in neutral matter. Therefore we impose that the total electric 
charge should vanish 


2 1 


0 . 


(36) 


The CSL condensates are color neutral such that no color chemical potentials 
are needed and no further constraints need to be obeyed. 

We calculate the CSL gaps for each flavor as a function of the quark chemical 
potential fx = {fXu + , 


^/-'^/(/^/(/^)) > ( 37 ) 

where the functional relation is taken from the (3- equilibrated and 

neutral normal quark matter equation of state. 

The results for the CSL gaps and the effective gaps as a function of the chem¬ 
ical potential fx are shown in Fig. 6 and in Fig. 7, respectively. From the 
parameterizations we listed in Table 1, we choose the Gaussian set because we 
consider it the most promising for stable hybrid star conhgurations due to the 
early onset of chiral and superconducting phase transitions in quark matter. 

From the hgures 6 and 7, we see that the two branches of the gap functions 
corresponding to the up and down quarks are put apart by the charge neu¬ 
trality condition (thick lines). The smallest gap, A^ runs from ^ 100 keV near 
the onset to ~ 500 keV at /i = 500 MeV while for the d quarks, A^ increases 
from ^ 380 keV to 1.4 MeV in the same range. 

On the other hand, the effective gaps Aj efr are of the order of magnitude of 
~ 10 keV, showing an approximate linear behaviour with /x. It remains to be 
investigated whether such small effective gaps could effectively suppress the 
direct Urea process in quark matter which is a requirement of compact star 
cooling phenomenology. 

In this respect, we found here that two facts produce a strong reduction of 
the CSL energy gaps. First, when we include smooth form factors in the ef¬ 
fective interactions, the values of the gaps decrease dramatically relative to 
the NJL case (an order of magnitude from ~ 100 keV for NJL to ~ 10 keV 
for Gaussian). Second, when neutrality constraints are considered, A„ and the 
effective gaps A/^eff are reduced further. 
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Fig. 6. CSL gaps for /3-equilibrated and neutral quark matter as a function of the 
quark chemical potential fi for the Gaussian set. 


However, as we have shown, there is a strong dependence of onr resnlts on 
the parameterization, and a more systematic investigation on the smoothness 
of the form factor could be helpful to decide whether these phases could be 
suitable for compact star applications. Moreover, this study should be seen as a 
preparatory step for subsequent investigations where, for example, a covariant 
generalization of the formalism for the inclusion of nonlocality effects [20,21] 
should be considered. 


6 Conclusion 


We have studied the effect of nonlocality on spin-1 condensates in two fla¬ 
vor quark matter: the 2SC-t-spin-l pairing of the blue quarks (2SCb) and the 
color spin locking (CSL) phase. We found that the size of these small gaps is 
very sensitive to the form of the regularization. The nonlocality has a strong 
impact on the low density region and we obtain an earlier onset for the su¬ 
perconducting phases. This might be crucial to stabilize quark matter cores 
in hybrid stars. 
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Fig. 7. Effective CSL gaps for /3-equilibrated and neutral quark matter as a function 
of /i for the Gaussian set. 


On the other hand, due to the flavor asymmetry, we And that the 2SCb pairing 
phase can be ruled out for compact stars applications. The CSL phase, in 
contrast, is flavor independent and therefore inert against the constraint of 
electric neutrality. For electrically neutral quark matter in beta equilibrium 
we obtain effective CSL gaps which are of the order of magnitude of 10 keV, 
which might help to suppress the direct Urea process, in accordance with 
recent results from compact star cooling phenomenology. Nevertheless, since 
our results are strongly dependent on the parameters used and on the form 
of the regularization, more systematic studies are needed in order to decide 
whether the CSL phase could be applied in the description of compact stars. 
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